Spin-orbit coupling induced polarization textures in trapped Fermi gases 
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We use the Bogoliubov-de Gennes (BdG) theory to study the interplay between the Zeeman field, spin-orbit 
coupling (SOC) and resultant skyrmion-like polarization textures in harmonically trapped two-dimensional 
Fermi gases. Our primary finding is that, in sharp contrast to the no-SOC case where only an easy-axis po- 
larization is possible beyond a threshold Zeeman field (depending on the interaction strength), any non-zero 
combination of the Zeeman field and Rashba SOC induces not only an easy-axis polarization everywhere in the 
system but also a spatially-modulated transverse one near the edges. Our secondary finding is that the ground 
states of Rashba spin-orbit coupled systems are revealed to have stronger Fulde-Ferrell type non-uniform su- 
perfluid characters (i.e. phase modulations) but not a Larkin-Ovchinnikov one (i.e. amplitude modulations) 
under the in-plane Zeeman field, for which the phases of the superfiuid order parameters modulate near the 
edges and most prominently along the transverse direction. The effects of the interaction, temperature, SOC 
anisotropy and Zeeman field anisotropy on the intricate polarization textures and superfiuid order parameter are 
also discussed. 

PACS numbers: 05.30.Fk, 03.75.Ss, 03.75.Hh 



I. INTRODUCTION 



The possibility of simulating non-Abelian artificial gauge 
fields with quantum Bose and Fermi gases in atomic sys- 
tems has become one of the forefront research directions in 
the atomic and molecular physics community QlJ-[5|], primar- 
ily due to its direct connection to the topological phases of 
matter that has extensively been studied in the condensed- 
matter community in recent years Oq-^O. In particular, the 
exciting possibility of the creation and observation of Majo- 
rana bound states in topological insulators, superconductors 
and superfluids is at the heart of topological quantum com- 
putation IU0I1 . These novel quasiparticles can be created at 
the boundaries (edges) of non-Abelian topological phases, and 
they allow for non-local storage of quantum information that 
is protected from local perturbations by the bulk gap. Moti- 
vated by these theoretical proposals, spin-orbit coupled Fermi 
gases have recently been created and detected near the quan- 
tum degeneracy limit by two groups 03101. While the Shanxi 
group in China studied the spin dynamics and momentum dis- 
tribution asymmetry in the equilibrium state as hallmarks of 
the spin-orbit coupling y|], the MIT group used a more di- 
rect approach and analyzed the energy-momentum dispersion, 
spin-orbit gap and spin composition of the quantum states [4]. 
Thus, assuming that sufficiently low temperatures are experi- 
mentally attainable in the near future, the physics of Majorana 
bound states can be studied in the clean and controllable envi- 
ronment uniquely offered by the atomic systems Illll[l2ll . 

Following the success of these initial experiments |Q]-|5||, 
there has been growing theoretical interest in studying the 
single-, two-, few-, and many-body properties of spin-orbit 
coupled Fermi gases. For instance, the stability and phase 
diagrams have been studied for finite and infinite (thermo- 
dynamic) systems as functions of the interaction strength, 
SOC strength, population imbalance, Zeeman fields, SOC 
anisotropy, Zeeman field anisotropy, temperature, etc. in one, 
two and three dimensions (see e.g. II13T - U8IO . The amount of 
knowledge gained in these recent works is overwhelming, and 



here we briefly quote the most recent ones that are concen- 
trated on the possibility of creation and observation of Fulde- 
Ferrell-Larkin-Ovchinnikov (FFLO) |HH HI type spatially- 
modulated non-uniform superfiuid phases under in-plane Zee- 
man field li 1 31 - 4 1 811 - We note that the FFLO-type phases in 
spin-orbit coupled Fermi gases are stabilized mainly by the 
asymmetry of the Fermi surfaces in momentum space, and this 
mechanism is in contrast with that of the condensed-matter 
ones where they are stabilized by the symmetric Zeeman mis- 
match in momentum space. 

In this paper, we study the interplay between the Zee- 
man field, SOC and resultant polarization textures in two- 
dimensional Fermi gases J2ll - l23ll . Our primary finding for 
harmonically trapped systems is that any non-zero combina- 
tion of the Zeeman field and Rashba SOC induces not only 
an easy-axis polarization everywhere in the system but also 
a spatially-modulated (ring-shaped in magnitude) transverse 
one near the edges. This is in sharp contrast with the trapped 
systems with no-SOC (and also with the thermodynamic sys- 
tems with SOC) where only an easy-axis polarization can be 
induced beyond a threshold Zeeman field. We show that the 
induced skyrmion-like [a |24J, |25[] polarization textures are 
finite-size edge effects, broadened by the harmonic trapping 
potential, and that their microscopic origin can be traced back 
to the counter-flow of spontaneous spin currents in the case of 
Rashba SOC. We also argue that the transverse polarization 
textures can be used to probe and characterize the topological 
phase transitions and the associated Majorana bound states in 
finite spin-orbit coupled Fermi gases. 

Furthermore, our secondary finding is that while the ground 
states of the spin-orbit coupled systems may have weak Fulde- 
Ferrell (FF) [19] type non-uniform superfiuid characters (i.e. 
phase modulations) but not a Larkin-Ovchinnikov (LO) Il20ll 
one (i.e. amplitude modulations) under the out-of -plane Zee- 
man field, the FF character of the superfluids is stronger for 
the Rashba SOC under in-plane Zeeman field. The FF-type 
phase oscillations are most prominent along the direction that 
is perpendicular to the Zeeman field. Since the superfiuid or- 



der parameters modulate only towards the edges (but not close 
to the center) of the system, where the densities of fermions 
are low and the magnitudes of the order parameters are small, 
it may be difficult to detect these modulations in atomic sys- 
tems at finite temperatures. We also comment on the effects 
of interaction, temperature, SOC anisotropy and Zeeman field 
anisotropy on the polarization textures and superfiuid order 
parameter. 

The rest of this paper is organized as follows. In Sec. |TTJ 
first we introduce the mean-field BdG Hamiltonian and then 
derive the self-consistency equations for the superfiuid order 
parameter, total number of fermions, and out-of- and in-plane 
spin polarizations. We numerically solve the resultant equa- 
tions and discuss the obtained results in Sec. [HI] Finally, the 
conclusions of this paper are briefly summarized in Sec. |IV] 



II. BOGOLIUBOV-DE GENNES THEORY 

The results mentioned above are obtained within the real- 
space BdG theory as we describe next. The mean-field Hamil- 
tonian for a two-dimensional Lx L square lattice can be com- 
pactly written as a 4L 2 x 4X 2 matrix-eigenvalue problem IJ26I1 . 
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where u a ni and v° xi are the components of the nth quasiparticle 
wave function at site i, a = {t, 4-} labels the z-component of 
the fermion spin, and e n > is the corresponding energy 
eigenvalue. The offsite hopping and onsite energy terms can 
be compactly written as 

+ (h x - is a h y ) 6 8crt - s<r ,] Sij, (2) 

where /x is the chemical potential, s-f- = — 84, = 1, dij is 
the Kronecker delta, h = (h x , h y ,h z ) is the Zeeman field, 
and Vi — Vbr 2 is the harmonic confining potential where the 
distance ?*j of site i is measured from the center of the lat- 
tice. The non-vanishing nearest-neighbor hopping elements 
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the positive x direction, and t"^ = icos/3 and £*? = 

t l ?l = it sin /3 for the positive y direction. Note that the 
hopping in the negative directions are simply the Hermitian 
conjugates, and also that the angles a and j3 determine, re- 
spectively, the relative strength between the spin-conserving 
particle hopping and spin-flipping SOC terms in the x and y 
directions. 

In this paper, we consider only the onsite interactions for 
which the off-diagonal couplings are A^ = AiSij diagonal 
in the site index, where the complex number Aj is the lo- 
cal mean-field superfiuid order parameter. Therefore, Eq. dD 
needs to be solved simultaneously with Aj = g(ctiC±i), 



where g > is the strength of the onsite interaction between 
t and I fermions, c* ai (c^) creates (annihilates) a a fermion 
at site i, and (• • • ) is a thermal average. In addition, we use 
/1 to fix the total number of fermions N — J2i n i> where < 
ni = X)ct( c L c '") — 2 gives the local fermion filling. Once 
the self-consistent solutions are obtained for the wave func- 
tions and the energy spectrum, it is a straightforward task to 
calculate any of the desired observables. For instance, we are 
interested in the local polarization vector pi = {pi X ,Piy,Piz), 
the components of which follow from the expectation values 
of the Pauli spin matrices, and are given by pi X = 2Re(cLc|,j), 

Piy 

the following averages 



= 2Im(cLc4.i) andp iz = ^ CT s G ^ ai c ai ). Thus, we need 
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(cliCai) = E [l<«! 2 /(e„) + K t \ 2 f(-e n ) 



(3) 
(4) 

(5) 



where f(x) = l/{e x / T + 1) is the Fermi-Dirac distribution 
function, T is the temperature and the Boltzmann constant 
ks is set to unity. We also define the total polarization com- 
ponents as P„ = J2iPiv/N , where v = {x, y, z}. Equa- 
tions ([T])-© correspond to the generalization of the BdG equa- 
tions to the case of spin-orbit coupled Fermi gases on optical 
lattices. 



III. NUMERICAL RESULTS 

Having established the BdG formalism, next we present 
our numerical solutions for the ground-state phases, which 
are performed on a 41a x 41a square lattice with TV = 150 
fermions in total, where a is the lattice spacing. We take 
Vq = O.Qlt as the strength of the trapping potential, and dis- 
cuss both the Rashba-type symmetric {a = /3) and asymmet- 
ric {a 7^ ft) SOC fields. Note that the experimentally more 
relevant equal Rashba-Dresselhaus (ERD) SOC [1-5] can be 
obtained by setting a = 0. The effects of higher fermion num- 
bers and finite temperature are also briefly mentioned towards 
the end of the paper. 

Before we discuss the aforementioned polarization textures, 
we make three important remarks. First, in the absence of a 
SOC, i.e. when a = f3 = 0, we know that a sufficiently strong 
Zeeman field h (the threshold of which depends on g) can po- 
larize the system along the easy-axis (h) direction. Second, 
in the absence of a Zeeman field, i.e. when h = (0, 0, 0), we 
also know that the system is trivially unpolarized no matter 
what the SOC is. Third, while any combination of Zeeman 
field (no matter how weak the field is) and SOC in a thermo- 
dynamic system may produce a uniform polarization along 
the easy-axis direction, it does not induce any polarization in 
the transverse direction, i.e. perpendicular to h. With these 
remarks in mind, next we show that any non-zero Zeeman 
field can induce intricate polarization textures near the edges 



of finite-size spin-orbit coupled Fermi gases under various cir- 
cumstances. 



A. Out-of-plane Zeeman field 

Let us first consider an out-of-plane Zeeman field h = 
(0, 0, h z ^ 0), which is perpendicular to our square lattice. 
The magnitudes and phases of typical ground-state order pa- 
rameters are illustrated for the Rashba and ERD-like SOCs in 
Figs. Q2a) and Hlb), respectively, for h z — 0.5i. The phases 
of the order parameters clearly show the d and C-2 symme- 
tries of the Hamiltonian for the Rashba and ERD-like SOCs, 
respectively. These figures also suggest that the ERD-like 
SOC has a stronger FF-type non-uniform superfiuid charac- 
ter where the phase of the order parameter has a much larger 
spatial modulation. Note that the phases have both angular 
and radial oscillations towards the edges of the system where 
the densities of the fermions are low and the magnitudes of 
the order parameters are small. Therefore, it may be diffi- 
cult to detect these modulations in atomic systems at finite T. 
However, since the spatial profiles of the magnitudes of the 
order parameters do not have any zeros (nodes), our results 
do not feature any LO-type non-uniform superfluidity. We 
emphasize that these effects become weaker and weaker with 
decreasing h z in such a way that all of the local phases vanish 
as h z -^ 0. 
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FIG. 1: (Color online) The color maps of the amplitudes (|Aj| - left 
column) and phases (tj>i - right column) of the order parameters Aj = 
| Ai |e * are shown for an out-of-plane Zeeman field h = (0, 0, 0.5i) 
at g — 3£ and T — 0, where a — j3 = 7r/4 in (a), and a = 7r/40 
and /3 = 7r/4 in (b). 

The corresponding ground-state polarization textures are il- 
lustrated in Fig. |2] where we show two-dimensional vector 
maps of the transverse polarizations (~Pi X ,~Pi y ) together 
with color maps of the easy-axis polarizations pi Z . Here, we 
set g = 3t but emphasize that setting it to does not lead 



to any significant change in the results. These figures again 
clearly show the C4 and C2 symmetries of the Hamiltonian for 
the Rashba and ERD-like SOCs, respectively. First of all, the 
not-so-interesting pi Z is finite everywhere in the trap with its 
maximum value at the center of the system in both figures, and 
it gradually decreases to zero towards the edges. The case of 
Rashba SOC is shown in Fig. [2 a), where we find that pi X ^ 
and pi y 7^ in general, except for the center of the trap. In the 
ERD-like case when a — > but j3 — 7r/4, we see in Fig.|2jb) 
that while p ix — )• everywhere in the system, p iy remains 
mostly unchanged. Therefore, in the ERD case when a = 0, 
a domain-wall is formed on the x axis where pi X = p- ly = 0, 
and such a limiting behavior can be extracted from Fig. EJb). 
Similarly, when /3 = but a ^ 0, a domain-wall forms on 
the y axis where pi X = pi y = (not shown). Thus, we con- 
clude that a non-zero out-of-plane Zeeman field no matter how 
small it is (not shown) induces spatially modulated transverse 
polarizations only along those directions where there is SOC. 
In most cases, the ratio of the transverse to the easy-axis po- 
larizations are around %5 — %10. However, we emphasize 
that while the total easy-axis polarizations are P z » 0.3 and 
P z « 0.36 in Figs. |2la) and[2{b), respectively, the total trans- 
verse polarizations vanish, i.e. P x = = P y , as one may 
expect. 

Note that the transverse polarizations change sign in space, 

and that their magnitudes \ p\ x + pf show ring-shaped struc- 
tures. In Fig.|2fa), in addition to the broad ring ranging mostly 
between la < r; < 12a, there is also a narrow one with a very 
weak peak around r, ps 14a. Similarly, in Fig. EJb), there 
are two incomplete rings especially along the y axis, and they 
both have comparable peaks around r, « {10.5a, 14.5a}. 
We also find that increasing the size of the square lattice 
pushes the ring-shaped structures further away from the cen- 
ter, and reducing the strength of the trapping potential inten- 
sifies them around a narrower region near the edges, eventu- 
ally both leaving no transverse polarization near the center. In 
the box-potential limit when Vq —} 0, we find that while the 
continuum- and edge-state contributions to the transverse po- 
larizations are competing with each other for low h z values, 
the latter contribution gets stronger with increasing h z , and 
eventually dominates beyond the h z threshold for the creation 
of zero-energy (Majorana) edge-bound states. These find- 
ings suggest that the transverse polarizations observed here 
are finite-size edge effects, broadened by the harmonic trap- 
ping potential, and also that the working mechanism is similar 
to the one that is responsible for the creation of edge-bound 
Majorana states. 

To understand the microscopic origin of these textures, next 
we employ the local-density approximation and analyze the 
single-particle excitation spectrum of the local system. The 
spectrum of the local Hamiltonian in momentum space in- 
volves two quasihole and two quasiparticle branches that are 
determined by 
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where £ki = ~2t[cosacos(k x a) + cos j3 cos(k y a)] — fit is 
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FIG. 2: (Color online) Two-dimensional vector maps of the trans- 
verse polarizations (— pi x , —pi y ) are shown for an out-of-plane Zee- 
man field h 5E (0, 0, 0.5t) at g = 3t and T = 0. Here, a = /3 = tt/4 
in (a), and a = 7r/40 and /3 = 7r/4 in (b). The largest arrows corre- 
spond approximately to 0.01, and the easy-axis polarizations pi Z are 
illustrated with color maps. 



the shifted kinetic energy and |sk| 2 = 4£ 2 [sin 2 a sm 2 (k x a) + 
sin 2 /3sin 2 (fcj,a)] is the SOC contribution. Here, the local 
chemical potential /ii = jX—Vi includes the trapping potential. 
We immediately see that the minus branches can become gap- 
less at some k-space points, i.e. £"£ j _ = 0, and therefore the 
location of zero-energy states are determined by the following 

I A,- 1 2 . While 



conditions: (i) |s ko | = 0md(ii)h z = J(i oi 

the Rashba SOC satisfies the former condition at four points 
k = {(0, 0); (0, tt); (tt, 0); (tt, tt)}, the ERD SOC satisfies it 
only at two points ko = {(fc x ,0); (k x ,ir)}. 

It is clear that the latter condition (ii) is easier to sat- 
isfy towards the edges of the system when |A,-| — > 0, and 
since the transverse polarizations are found to be very sim- 
ilar for g = 3t and g = 0, we may set |A,-| = for our 
purpose. When this is the case, the condition (ii) becomes 
fii = ±h z ± 2t(cosa ± cos/3) for the Rashba SOC, and 



Hi = ±h z — 2t[cos(k x a) ± cos/3] for the ERD SOC, where 
all ± combinations are possible and | cos^a)! < 1. For the 
parameters of Figs. 12 a) and|2jb), where /i ra — 1.76t, these 
conditions are satisfied at two distances r.i k. {7.5a, 12.5a} 
and Ti « {10.7a, 14.6a}, respectively, which are very close 
to our numerical results given above. Thus, we conclude that 
the microscopic origin of the transverse polarizations can be 
traced back to the changes in the momentum-space topology 
of the single-particle excitation spectrum of the local system. 
In the case of Rashba SOC, we can also interpret 
these textures as a direct consequence of counter-flow of 
spontaneously-induced spin currents [26]. This is because the 
Rashba SOC gives rise to an effective momentum-dependent 
in-plane magnetic field in the direction that is perpendicular to 
the in-plane momentum. Since the induced spin currents are 
circulating along the trap edges, i.e. the in-plane momentum 
is in the azimuthal direction, the induced in-plane spin texture 
is in the radial direction. The relative contribution between 
the radially-outward and -inward helicity bands depends on 
the local chemical potential fa, and this competition produces 
the spatial structure of the spin textures such as the one illus- 
trated in Fig. 12 a). Note that the time-reversal symmetry of 
the spins must be broken via e.g. the Zeeman field in order to 
have a non-zero polarization in any particular direction. 



B. In-plane Zeeman field 

Next, we consider an in-plane Zeeman field h = (0, h y j^ 
0, 0), which lies in the y direction parallel to our square lat- 
tice. The magnitudes and phases of typical ground-state order 
parameters are illustrated for the Rashba and ERD-like SOCs 
in Figs. O a) and|3jb), respectively, for h y = 0.5t. In sharp 
contrast to the out-of-plane Zeeman case discussed above, this 
comparison clearly shows that the Rashba SOC has a much 
stronger FF-type non-uniform superfluid character in the in- 
plane Zeeman case, without again featuring any LO-type or- 
der parameter node. Note again that the phases have both an- 
gular and radial oscillations towards the edges of the system, 
and the FF-type oscillations are most prominent along the x 
direction, i.e. perpendicular to the direction of the Zeeman 
field. We emphasize that these effects become weaker and 
weaker with decreasing h y in such a way that all of the local 
phases vanish as h y — > 0. 

The corresponding ground-state polarization textures are il- 
lustrated in Fig. |U where we show two-dimensional vector 
maps of the transverse polarizations (pi x , pi z ) together with 
color maps of the easy-axis polarizations pi y . Here, we set 
g = 3t but setting g to again does not lead to any signifi- 
cant change in the results. First of all, the not-so-interesting 
Piy is finite everywhere in the trap with its maximum value 
near the center in both figures, and it gradually decreases to 
zero towards the edges. The case of Rashba SOC is shown in 
Fig. Ufa), where we find that pt x ^ and pi z ^ in general, 
except for a domain-wall on the x axis where pi x = p. LZ = 0. 
In the ERD-like SOC when a -»■ but /3 ^ 0, we see in 
Fig. Ub) that while pt x -> everywhere in the system, pi Z 
remains mostly unchanged. Therefore, similar to the Rashba 
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FIG. 3: (Color online) The color maps of the amplitudes (|Aj| - left 
column) and phases (4>i - right column) of the order parameters A; = 
lAile 1 ^* are shown for an in-plane Zeeman field h = (0, 0.5i, 0) at 
g — 3t and T — 0, where a — fi — n/4 in (a), and a = 7r/40 and 
j3 = tt/4 in (b). 



case, the ERD case also has a domain-wall that is formed on 
the x axis where pi X = pi Z = 0, and such a limiting be- 
havior can be extracted from Fig. HJb). On the other hand, 
when j3 = but a ^ 0, there is not any transverse polariza- 
tion in the entire system, i.e. pi X = pi Z = for every i (not 
shown). Thus, we conclude that, when a non-zero in-plane 
Zeeman field (no matter how small it is) is not perpendicular 
to the direction of the SOC, a spatially modulated polarization 
is induced in the transverse direction. In most cases, the ra- 
tio of the transverse to the easy-axis polarizations are around 
%5 — %15. However, we emphasize that while the total easy- 
axis polarizations are P y w 0.37 and P y ss 0.36 in Figs. Ufa) 
andQb), respectively, the total transverse polarizations again 
vanish, i.e. P x = = P z . 

Similar to the out-of-plane Zeeman case, we note that the 
magnitudes of the transverse polarizations \Jp 2 ix + p 2 z show 
ring-shaped structures in both Figs. |4ja) and|4jb), where the 
peaks occur, respectively, at r^ rs {10.5a, 13.5a} and r, ~ 
{10.5a, 14.5a} away from the center of the trap especially 
along the y axis. The microscopic origin of these structures 
can again be traced back to the changes in the momentum- 
space topology of the single-particle excitation spectrum. 



C. Topological phase transitions 

As elaborated above, although gap closings occur only at 
a few ko points in the lowest quasiparticle and highest quasi- 
hole bands, these gapless excitation points are sufficient to 
induce intricate polarization textures in real space, at and 
around the boundary between phases with locally different 
momentum-space topology. We emphasize that since the 
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FIG. 4: (Color online) Two-dimensional vector maps of the trans- 
verse polarization (j3ix,Piz) are shown for an in-plane Zeeman field 
h = (0,0.5t,0) at g = 3t. Here, a = /3 = tt/4 in (a), and 
a — 7r/40 and /? = n/4 in (b). The largest arrows correspond ap- 
proximately to 0.014, and the easy-axis polarization pi V is illustrated 
with a color map. 



symmetry of the order parameters of the local phases have 
the same s-wave symmetry across the boundary, the transi- 
tion is topological. Such topological changes are known as 
Lifshitz-type transition in the literature, and they are exten- 
sively discussed in the context of nodal, e.g. p-wave, super- 
fluids and superconductors [6]. In thermodynamic systems, 
while the primary signatures of Lifshitz transitions are seen in 
the momentum distribution and single-particle spectral func- 
tion, some thermodynamic quantities, e.g. atomic compress- 
ibility and spin susceptibility, also show anomalies at the tran- 
sition boundary. 

It is also worth noting that both spinless p x ± ip y (chiral) 
superfluids [6] and spin-orbit coupled Fermi gases under Zee- 
man field llll 11211 can host Majorana bound states near the 
edges of the system. More specifically, these bound states 
can only exist at the phase boundary between a topologically 
non-trivial and a trivial phase, the classification of which is 



based on the value of the topological charges, i.e. Chem num- 
bers Da]. In contrast to the spinless chiral superfluids, our nu- 
merical results suggest that the induced transverse polariza- 
tions can be used as a probe to characterize the topological 
phase transitions and the associated Majorana bound states in 
finite spin-orbit coupled Fermi gases. This is in accordance 
with a recent work on one-dimensional quantum wires with 
strong Rashba and Dresselhaus SOC, where it is shown that 
the Majorana polarization can be used as an order parameter to 
characterize the topological transition between the trivial sys- 
tem and the system exhibiting Majorana bound modes 112711 . 

Having discussed the low-filling Fermi gases at zero tem- 
perature, next we briefly comment on the effects of finite T 
and high fillings. First, although these topological transitions 
are quantum in their nature, signatures of them can still be 
observed at finite T, where the observables are smeared out 
due to thermal effects. In particular, we find for the param- 
eters of Figs. [2] and [4] that the maximum magnitudes of the 
transverse polarizations reduce, respectively, to 50% and 10% 
at T = O.lt and T = 0.2t. Second, due to the particle-hole 
symmetry of the parent Hamiltonian around half filling, in ad- 
dition to the ring-shaped transverse polarizations induced near 
the edges of the system, additional ring-shaped structures are 
further induced near the center of the trap when the center is 
close to a band insulator. Therefore, the transverse polariza- 
tions show multiple ring-shaped structures in high-filling lat- 
tice systems. Having discussed the numerical results, next we 
conclude the paper with a brief summary of our main findings. 



IV. CONCLUSIONS 



only an easy-axis polarization is possible beyond a threshold 
Zeeman field, we showed that any non-zero combination of 
the Zeeman field and Rashba SOC induces not only an easy- 
axis polarization everywhere in the system but also a spatially- 
modulated (ring-shaped in magnitude) transverse one near the 
edges. We found that the induced skyrmion-like polarization 
textures are finite-size edge effects, broadened by the har- 
monic trapping potential, and that their microscopic origin 
can be traced back to the counter-flow of spontaneous spin 
currents in the case of Rashba SOC. We also argued that the 
transverse polarizations can be used to probe and characterize 
the topological phase transitions and the associated Majorana 
bound states in finite spin-orbit coupled Fermi gases. 

Second, we showed that while the ground states of the 
spin-orbit coupled systems in general have weak FF-type non- 
uniform superfluid characters (i.e. phase modulations) but 
not an LO one (i.e. amplitude modulations) under the out- 
of -plane Zeeman field, the FF character of the superfluids is 
stronger for the Rashba SOC under in-plane Zeeman field. 
The FF-type phase oscillations are also most prominent along 
the direction that is perpendicular to the Zeeman field. Since 
the superfluid order parameters modulate only towards the 
edges (but not close to the center of the trap), where the den- 
sities of fermions are low and the magnitudes of the order pa- 
rameters are small, it may be difficult to detect these modula- 
tions in atomic systems at finite T. Furthermore, we discussed 
the effects of interaction, temperature, SOC anisotropy and 
Zeeman field anisotropy on the aforementioned polarization 
textures and superfluid order parameter, and briefly discussed 
the possibility of observing these effects in atomic systems. 



In this paper, we studied the interplay between the Zeeman 
field, SOC and resultant polarization textures in harmonically 
trapped two-dimensional Fermi gases. The trapping poten- 
tial is taken into account via the real-space mean-field BdG 
theory, and two of our main findings can be summarized as 
follows. First, in sharp contrast to the no-SOC case where 
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